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The paradigm of canonical Hamiltonian perturbation theory Q involves the transfor- 
mation of an exact Hamiltonian H^, which depends continuously on a perturbation param- 
eter e, into a reference Hamiltonian Hq for which the Hamilton equations dtz = {z, Hq} 
have a known solution (unless otherwise noted z denotes canonical phase-space coordinates 
and { , } denotes the canonical Poisson bracket). According to the Lie-transform approach 
to Hamiltonian perturbation theory the transformation — »• Hq is induced by a re- 
versible (time-dependent) phase-space transformation from the old phase-space coordinates 
z to the new phase-space coordinates z(e) = z, where is an operator defined in terms 
of a generating scalar field as 



T, = exp (^J^ {S„, } da 



For a time-dependent phase-space transformation z z(e) generated by the scalar field 
S^{z,t), the transformation from the old Hamiltonian to the new Hamiltonian Hq is 
expressed as 0, ^ 

i/o(z(6),t) = H,iz,t) - r ^^^^rfa. (1) 

Jo ot 

The transformations z —>■ z and if^ ^ iJo are therefore completely determined by the 
perturbed Hamiltonian H^{z,t) and the phase-space generating function S^{z,t). Since 
Hq is independent of e by construction (i.e., 9e-?fo = 0), the e-derivative of both sides 
in (HI) yields a dynamical evolution (henceforth known as the Lie-transform perturbation 
equation) 

dtS + {S, H} = d,H, (2) 

where the parametric e-dependence is included with the time dependence. The e-perturbed 
Hamilton equations, on the other hand, are now expressed as 



dz''{t,e) 
di 

and 

dz''{t,e) 



H{z;t,e)}, (3) 



{^", S{z;t,e)}. (4) 



de 

Hence, whereas the Hamiltonian H is the generating function for the infinitesimal canonical 
transformations associated with the dynamical evolution (^ of the Hamiltonian system 
(henceforth referred to as the t-dynamics), the generating function 5* is the Hamiltonian 
for the perturbation evolution (^ of the Hamiltonian system (henceforth referred to as 
the e-dynamics). We note that the order with which the Hamiltonian system is evolved 
and perturbed should be immaterial, i.e., the same dynamical state z{t,e) = z^(t) can 
be reached by either evolving the unperturbed system first [zo(t = 0) — > zo(t)] and then 
perturbing it [zo(t) ^^(t)] or perturbing the system first [zo(t = 0) z^(t = 0)] and 
then evolving it [z^(t = 0) — * ^^(t)]. 
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The purpose of this Letter is to present the variational formulation for the Lie-transform 
perturbation equation (|D and the multi-Hamilton equations (j^) and @. In particular, we 
show that the Lie-transform perturbation equation (0) is a direct consequence of path 
independence in our variational formulation. First, we introduce an extended phase-space 
Lagrangian T on the (t, e)-plane defined as 

r(Z(t,e);t,e) = P{t,e) • dQ{t,e) - H{Z{t,e);t,e) dt - S{Z{t,e);t,e) de, (5) 

where 

Z": (t,e) ^ Z"(t,e) = (Q(t,e),P(t,e)) (6) 
denotes a generic mapping from the {t, e)-plane to the 2A^-dimensional phase space, with 

ot oe 

(In what follows, the uppercase Z denotes a generic mapping whereas the lowercase z 
denotes a Hamiltonian orbit in phase space.) The one-form (H) is said to be extended in 
the sense that the e-dynamics Hamiltonian term — S de has been added to the standard 
phase-space Lagrangian P -c/Q — H dt. Next, we define the action integral 

Acm ^ I r(Z(t,e);t,e), (7) 

where C denotes an arbitrary path between two (distinct) points on the two-dimensional 
(t, e)-plane. For a fixed path C, the action integral ^c[Z] is a functional of the mapping 
(^. Holding C fixed, we first consider the variational principle ^^^[z] = corresponding 
to an arbitrary variation 5Z = Z — z (which is assumed to vanish at the end points of C). 
Using Eqs. (^ and (0), we thus find 

5AM = I 5Z-(^u^,dz^ - - ^^') + L ^(P-^Q)' (8) 

where the second path integral vanishes since we assumed that 5Z vanishes at the end 
points of C; here, ujap denotes the components of the canonical Lagrange tensor (i.e., 
dpi A dq^ = \ uja(3 dz" A dz^). The multi-Hamilton equations (|) and (0) are automatically 
recovered in Euler-Lagrange form as dz^ = daH dt + d^S de from the variational 
principle ^^ciz] = for arbitrary variations SZ and an arbitrary path C. 

In the standard variational principle for single-particle Hamiltonian dynamics [Q, the 
action functional involves a time integration of the phase-space Lagrangian P • dQ — H dt 
from an initial time ti to a final time the initial and final times play no role and the 
variational principle yields the standard Hamilton equations. The variational principle 
based on the action functional (|^ presents us with a new problem: determining how the 
freedom of choice in selecting the path C on the {t, e)-plane is expressed mathematically. 
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To resolve this problem, we first consider two different paths C and C, both having the 
same end points, and we choose the mapping (Q) to be a multi-Hamiltonian orbit z(t,e), 
i.e., a solution of Eqs. and (^). Next, according to Stokes' theorem [Q, the difference 
between ^c[z] and ^c'[^] is evaluated as 



Ac[z] - Ac'[z] = f T = dT, (9) 

Jc~C' JV 

where T> denotes the area in the (t, e)-plane enclosed by the two paths C and C (i.e., 
C — C = dV denotes the contour of V). Using Eqs. (0)-(|^), the two-form dT appearing in 
Eq. (H) is 

rfr(t,e;z(t,e)) = de A dt ( ^ - ^ - {S, H}] , (10) 

\ dt de J 

where the operators d/dt and d/de are defined as 

d/dt = dt + {, H} 

d/de = d, + {, S} 



(11) 



From Eq. (^), we see that, for a fixed Hamiltonian orbit z{t, e), the condition for the path- 
independence of the action integral ^c[z] (i.e., ^c[z] = ^c"[^]) is dT = 0. From Eq. ([T0|), 
we see that this condition holds provided H and S satisfy the following constraint equation: 

which is exactly the Lie-transform perturbation equation (§) appearing in Lie-transform 
Hamiltonian perturbation theory. We can also verify by using the operators defined in (pJ]) 
and the Jacobi identity for the Poisson bracket that the constraint ([l2|) ensures that the 
flows associated with t-dynamics and e-dynamics commute: 





= ^ Pf^ - = (13) 



d d 
dt^ de 

for all g{z;t,e). As expected the interpretation of the commutation relation (pTSj) is that 
the order with which the Hamiltonian system is evolved (under the t-dynamics) and is 
perturbed (under the e-dynamics) is indeed immaterial. 

So far we have presented the variational formulation of Lie-transform canonical Hamil- 
tonian perturbation theory. Most recent applications of Hamiltonian perturbation theory 
in plasma physics are carried out using non-canonical phase-space coordinates [^, |^ . Non- 
canonical Lie-transform Hamiltonian perturbation theory possesses a variational formula- 
tion expressed in terms of the extended non-canonical phase-space Lagrangian 



e ^ , , 
mv + - A(x; t, e) 
c 



' (ix — H{yi,\;t,e) dt — S{'x,v,t,e) de, (14) 
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where the non-canonical coordinates (x, v) denote the particle position and its velocity 
while A(x; t, e) denotes the perturbed magnetic vector potential. For a fixed path C in the 
{t, e)-plane, the variational principle 6{Jq F) = yields the Euler-Lagrange equations 

9H , dS , 
md'^ = — — at + — — ae 
av av 

mdv = - VH H —]dt- V5 H — ] de + - rfx x B, 

y c ot J \ c oe J c 

which can be written in multi-Hamilton form as 

dtz^ = {z", H}n, + (e/c)9iA.{x, ^ 

(15) 

d,z^ = {z-, 5}nc + (e/c) d,A ■ {x, ^ 
where { , }nc denotes the non-canonical Poisson bracket 

= 1 fva — - — \/b] + (— X —] 

The condition dT = for path independence, on the other hand, yields the Lie-transform 
non-canonical perturbation equation 

— - — = {H S} + —(—.— - —.—] (16) 
dt de ' "'^ mc \ dt dv de dv J 

This non-canonical Hamiltonian perturbation equation plays a prominent role in the deriva- 
tion of the reduced nonlinear gyrokinetic equations describing the perturbed Hamiltonian 
dynamics of charged particles under the infiuence of low-frequency electromagnetic fiuctu- 
ations in a magnetized plasma p. A linearized version of (16) also plays a prominent role 
in the free-energy method developed by Morrison and Pfirsch 0, H, ^ to investigate the 
linear stability of various plasma equilibria. 

Finally, we note that the variational formulation for Lie-transform Hamiltonian pertur- 
bation theory can be generalized to include several perturbation parameters. Consider the 
{k + l)-component perturbation vector e = (e°, e^, e^), where each perturbation parame- 
ter is paired with a Hamiltonian Sa (with e° = t and 5*0 = H). Hence, the multi-Hamilton 
equations (|^) and (^ are now replaced by the set of multi-Hamilton equations: 

= {^^ 5„(z;e)}, (17) 

where a = 0, 1, k. We now define the action functional Ac['^] = Jc r(Z(e); e) in terms 
of a path C in the {k + l)-dimensional e-space and the one-form: 

F(Z(e);e) = P(e) • dQ(e) - 5,(Z(e); e) rfe^ (18) 
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where Z : e i— Z(e) denotes a generic mapping from e-space to phase space. For a fixed 
(but arbitrary) path C, we recover the muhi-Hamilton equations ([T7|) from the variational 
principle = 0. In addition, for a fixed multi-Hamiltonian orbit z(e) = (q(e), p(e), 

the {k + 1) equations corresponding to the condition for the path-independence of ^c;[z] 
(i.e., dr = 0) are 

daSt - dtSa = {Sa, 5J. (19) 

Note that this set of equations is invariant under the gauge transformation Sa ^ Sa — 
dxl de'^, where x(q; e) is an arbitrary scalar field. Gauge invariance of the multi-Hamiltonian 
dynamics requires that the phase-space Lagrangian F transform as F —>■ T + dx and 
p ^ p + dx/dq. The gauge properties of the generating function S are further eluci- 
dated in Ref. ITO|. 

The work presented here introduced a variational formulation for Lie-transform Hamil- 
tonian perturbation theory. This work follows recent developments in the variational for- 
mulation of exact and reduced Vlasov-Maxwell equations |jTl], 0, |12|, |13[ as well as other 



important dynamical field equations in plasma physics (e.g., kinetic-MHD equations |]TU 



and drift- wave equations [14 
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